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Abstract 

We use an analytic perturbation expansion for the two-band system of tight 
binding electrons to discuss Bloch oscillations and Zener tunneling within this 
model. We make comparison with recent numerical results and predict ana- 
lytically the frequency of radiation expected from Zener tunneling, including 

its disappearance, as a function of the system parameters. 
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One of the earliest predictions of the theory of electrons in periodic crystals is of Bloch 
oscillations — time periodic motion at frequency ujb = Fd/h under application of a uniform 
constant field F in a crystal potential of lattice period d, associated with Bragg reflec- 
tion at Brillouin Zone boundaries. But in ordinary atomic crystals, subjected to accessible 
applied (electric) fields, scattering typically disrupts the coherence of electronic motion in 
times very short compared to the Bloch period Tb = 2tt/ujb, preventing observation of the 
phenomenon. However, the lattice constant d of semiconductor superlattices, typically two 
orders of magnitude larger, correspondingly reduces Tb by the same factor. Esaki and TsJ 
recognized 25 years ago, at the dawn of this new technology, that this greatly improved the 
prospects of experimental observation of Bloch oscillations and related phenomena. More- 
over, the characteristic frequencies lie in the interesting far infrared region. Not surprisingly, 
then, there has been a great deal of subsequent activity, both experimental! and theoretical!, 
in the study of high quality semiconductor superlattices under the influence of static and 
time periodic electric fields. Transport and both linear and nonlinear optical response have 
been of special interest. In particular, Bloch oscillations have been shown! to survive the 
addition of Coulomb interactions to the simple theory, and far infrared radiation from such 
oscillations has now@ been observed. 

The one- dimensional periodicity of a quantum well superlattice gives as the eigenstates 
for independent (noninteracting) electrons a set of "minibands" associated with motion in 
the growth direction, perpendicular to the wells. For obvious reasons of simplicity, much of 
the theoretical study has been limited to consideration of a single isolated miniband. But, 
as emphasized particularly in a recent Letter! by Rotvig, Jauho and Smith (RJS), there 
are essential interesting phenomena — notably Zener tunneling — introduced by additional 
bands in the presence of a static electric field. These authors reported a numerical study of 
the density matrix of a two-band model, exhibiting a number of interesting types of periodic 
behavior of the electrons. It is the purpose of the present paper, (i) to make use of analytic 
results in the form of an exact perturbation expansion to explore the corresponding behavior 
throughout the space of the parameters characterizing the two-band superlattice, and (ii) to 



make use of our analytically predicted behavior to compare with selected numerical results in 
order to test the convergence of the perturbation series so as to define its range of validity for 
application to other problems and properties. The structure of the series suggests (correctly, 
we will show) that the convergence is much more rapid than imposed by the obvious limits 
that can be set analytically. 

The two-band model is, of course, the simplest extension beyond the single band picture, 
and it does introduce those features characteristic of interband communication. But is it a 
reasonable approximation to any physical situation to focus on a single pair of minibands 
while neglecting the higher bands which inevitably also exist? As we have pointed out 
beforejll there is at least one realistic situation where this should capture the dominant 
physical behavior. A superlattice can be grown dimerized, with the unit cell consisting, e.g., 
of a pair of quantum wells, separated from the adjacent pair by a larger (wider) barrier than 
that which separates the two wells of the basis pair. Then the lowest level of an individual 
well is split within the "molecular" pair. These two levels form a pair of bands in the 
superlattice which are well separated energetically from those arising from higher levels in 
the well, and it is then reasonable to neglect all except that lowest pair of bands. 



We will use the standard two-band tight-binding Hamiltonian in a uniform static electric 
field E: 



Here the subscripts label the lattice sites and the lower and upper minibands are designated 
by symbols a and b, respectively. The first two terms describe the site energies of the Wannier 
states in the presence of the electric field, and W a> b are the widths of the isolated (E = 0) 
minibands induced by nearest neighbor hopping: e a ' b {k) = A a ^ =f (W a ,b/2) cos(kd). The last 
term is the on-site electric dipole coupling between minibands; eR is the corresponding dipole 
moment. This Hamiltonian, first used! by Fukuyama, Bari and Fogedby in connection with 



H = (A a + neEd)a ] n a n + (A 6 + neEd)b ] n b, 



-(W a /4)(ai +1 a n + h.c.) + (W b /4)(b[ +1 b n + h.c.) 




(1) 
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this type of system, does neglect Coulomb interactions and electric dipole elements between 
Wannier states on different sites, but it contains the essential physics for the problem. Note 
that the hopping parameters W a ,b are written here with opposite signs, so that with both 
parameters positive the band structure at E = is of the standard nearly free electron 
character, with direct band gaps at the zone boundary. But the calculation to follow is valid 
for arbitrary signs of the parameters. 

The application of a static electric field localizes all electronic energy eigenstates, with 
a localization length which decreases with increasing electric field strength E. The energy 
spectrum associated with the Hamiltonian ([I]) above becomes two interleaved Stark ladders: 



where u>b = eEd is the Bloch frequency defined above (we have taken units with Planck's 
constant equal to unity), m is an arbitrary integer, and the Stark ladder offsets ilx,2 depend 
on the field amplitude E. The result so far can be established in several waydli. A partic- 
ularly satisfying proof based on the symmetries evident in a suitably chosen gauge is given 
in the Appendix. But one of us has established,! in addition, an analytic solution which 
includes an explicit expression for the offset between the two ladders Q = Q 2 ~ ^i- As noted 
by RJS, the Zener tunneling behavior undergoes dramatic changes when Q approaches zero. 
We first analyze that interplay between ordinary Bloch oscillations and Zener tunneling 
with the use of those analytic results, to understand how the behavior is affected by the 
parameters of the superlattice. One of the most attractive features of the semiconductor 
systems is, after all, that those parameters are to a considerable extent under the control of 
the fabricator. 

Formally, the analytic results are in the form of a pertubation series in the ratio R/d 
of the interband dipole coupling to the Bloch frequency. It is therefore natural to compare 
the exact energy eigenvalues to those of the decoupled bands (for the Hamiltonian (TO) with 
R = 0). These are of the same form as (fj): 




Q li2 + muj B , 



(2) 



€q' (n) = A atb + nu B , 



(3) 
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but the offsets A a> b are now the bare energy parameters of the Hamiltonian, independent of 
the applied field E. These energies are plotted in Fig. 1, as in RJS, in the form (e^ — A a )/u B , 
which equals n (for p = a) or (A a — Aj^/ub + n (for p = b), as a function of (A a — A^/Vb, 
giving a set of horizontal lines at the integers and a set of parallel lines of unit slope for 
bands a and b, respectively. Away from the crossing points of this unperturbed spectrum, 
A;, — A a = nu>B, where the two Stark ladders are degenerate, the exact full energies are 
expected (and found) to show only small corrections due to finite interband coupling R for 
reasonable values of R/d < 1. Similarly, the eigenstates are largely those of individual bands 
a and b. Then significant Zener tunneling is to be found only near the crossing points, and 
we will concentrate on those special values of the field, 

lo b = eEd = (A 6 - A a )/£, (4) 

where I is an integer. At these fields the eigenstates will be mixtures of comparable weights 
of states in the two bands. If the system is prepared initially in one of the bands, Zener 
tunneling will exhibit oscillations of the weight between the two bands with a dominant 
frequency of the Stark ladder offset Q. We will show, moreover, that there are special 
circumstances where this offset vanishes and the tunneling disappears. 

The results of Ref . |] can be expressed for the resonance condition of Eq. @ in the form 

n = -^cos- 1 Y,(R/d) 2m u {2m) {e,(w b + W a )/u B ], (5) 

71 m=0 

where we have explicitly indicated the dependence of the terms [/( 2m ) of the perturbation 
series on the parameters of the system. The contribution f/( 2m ) (other than the first: = 
1) is a sum of 2 2m terms, each of which is an ordered integral of the form 

/ dki / dk 2 --- dk 2m F 1 (k 1 )F 2 (k 2 ) ■ ■ ■F 2m (k 2m ), (6) 

Jo Jo Jo 

where each function F n (k) is of one of the two forms: 

C(k) = cos(a sin k — 2£k) 

S(k) = sin(asmk-2£k). (7) 
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Here a = (Wj, + W a )/ (2ub), and the integer £ is again the resonance parameter of Eq. (f|). 
The formal derivation of the perturbation series shows that |£/( 2m )| < (2vr) 2m /(2 m)\, so 
the series is convergent for all values of the ratio R/d. On the other hand, the number of 
terms and their algebraic complexity increase very rapidly with index (2m), so the result 
is of practical use only when it converges so rapidly that one or two terms suffice. In that 
regard the simple limit just invoked, which puts an upper bound on the term of order n 
of (27r R/d) n /n\, is not very helpful for typical values of 2irR/d, which are of order unity. 
But the highly oscillatory character of the integrands suggests that the series converges, in 
fact, much more rapidly than this. Our results below verify that this is indeed the case. We 
will keep only the first nontrivial correction in R/d, namely the m = 1 term, for which the 
integrals can readily be expressed in closed form: 

=2[nJ e (a)} 2 , (8) 

where Jg(a) is the Bessel function of order i. Already we see the substantial reductions 
from the oscillatory integrands from the crude limit above, though we have been unable 
to establish rigorous lower analytic limits than this. Near x — 1, cos -1 x has a square 
root singularity; for (R/d) 2 U^ C 1 we have as the offset between closest steps on the two 
ladders, 

nd 

= 2 Ma) (9) 

We note, in particular, that this drops off rapidly with increasing I (with decreasing field 
E for fixed (A& — A a ), as defined in (Q)), or as we move to the right in the level crossings 
shown in Fig. 1. This can be observed in the numerical results of Fig. 1 of RJS. The 
parameters used there are: A b — A a = 20meV, R/d = — 16/97T 2 = —0.18, and a = 0.43£, 
giving for the first three gaps (£ = 1,2,3) the values Q = 1.51, 0.306, 0.030 meV. For 
more direct comparison with the figure we reduce this by the Bloch frequency in each case: 
Q/ujb = 0.076,0.031,0.015, which appears to be in good agreement with those numerical 
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results. 

How accurate, in general, is the theory cut off at this lowest nontrivial order? To answer 
that, as well as to understand more easily the nature of the Zener tunneling, it is useful 
to rederive the result (|9[) without the full machinery of the exact solution, but rather using 
standard degenerate perturbation theory. For the independent band problem (R = 0) the 
Wannier-Stark eigenstates for band a are given explicitly as 

C = X>n-m (^)4J0>, (10) 

with energy given by (|3[). The wave function for band b is, of course, of precisely the 
same form. The resonance condition (|) then gives as degenerate the states </>^ and (j) b n+ ^ for 
each integer n. The degeneracy is lifted by the interband dipole matrix elements, which are 
given explicitly between degenerate unperturbed states as 

= /bE^(^)^(^) = -jUBrt*). (11) 

Then the degenerate states are split by twice this amount in lowest order, in agreement with 
(^), and the corresponding eigenstates are the symmetric and antisymmetric combinations, 

rrt a = (€±€ + i)/V2. (12) 

Now the nature of the Zener tunneling is transparent. Let us prepare the system in the 
(lower) band a before turning on the field at time t = 0. Then the subsequent time devel- 
opment is given by 

n 

where the symbols, other than the arbitrary constant coefficients A n which determine the 
initial state, have all been defined above. Clearly, the state becomes completely made of 
contributions from band b after time t = tt/Q and returns to band a in twice that time - 
i.e., with a full cycle of frequency Q, as seen in the numerical results of RJS. In particular, 



at the zeroes of the Bessel function, Ji(a) = 0, the splitting of the ladders Q also vanishes, 
and Zener tunneling also disappears. Although this is tunneling between bands rather than 
between the spatially separated sides of a standard double well, the physical manifestations 
are not dissimilar. Corresponding degenerate states from the two bands are shifted by 
I lattice sites, and there will therefore be an ac current, or oscillating dipole associated 
with the Zener tunneling, and corresponding radiation at the frequency Q. Therefore, the 
prediction of the theory to this order is that such radiation will be observed at a frequency 
Q oc Je(a), except for the special points where the Bessel function vanishes, at which the 
radiation will disappear, as well. These special points for various resonance ratios I (electric 
field given by are given by 

£(W a + W b ) 

a =2(A^ = X - (14) 

where Xt> is a zero of Jt{x). 

The full theory predicts similar phenomena, but then the splitting of the Stark ladders 
fl needs to be calculated to all orders. If the corrections are significant, then the predictions 
of (0) and (|TJ) will be at least quantitatively incorrect. 
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APPENDIX A: PROOF OF INTERLEAVED STARK LADDER SPECTRUM 

The translational symmetry of the Hamiltonian (JJ) is broken by the scalar potential 
terms proportional to the intraband part of the position operator, 
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x° 



dY,n{ala n + b%). (Al) 

n 

It is convenient to restore that translational symmetry by a suitable gauge transformation 
(different from the standard vector potential gauge, because only the intraband part of the 
position operator is involved): 

which obeys the time-dependent Schrodinger equation (id t — Ti')^' = 0, with 

w = XX 

k 

H' k = [A„ - (W„/2) cos(kd - uj B t)]a\k)a(k) 



+ 



A b + (W 6 /2) cos(kd - u B t)] b\k)b{k) 



+eER[a^(k)b(k) + tf(k)a(k)} . (A3) 
Here a(k), the spatial Fourier transform of a n , 

a ( fc ) = ^Ee mfc V, (A4) 

destroys a plane wave state of wave vector k in band a, and the hermitian conjugate operator 
a\k) creates that state. The operators b(k) and b\k) are defined analogously for the other 
band. The dynamical equation thus separates into a set of independent two-level problems, 
one for each wave vector k, described by Ti' k . Moreover, Tif k (t + TB) = T-L' k (t), a discrete time 
translational symmetry which, by the Floquet theorem, implies solutions of the form 

$fc(*) = u k(t) exp(-ie k t), (A5) 

with u k (t) = u k {t + T B ) periodic in time, and the "quasienergy" e k defined only modulo uj b . 
Further, the form ( |A3| ) of TC' k makes it clear that e k and u k are, in fact, independent of k: 
the parameter k merely shifts the origin of time, t' — t — kd/u B in each term. Therefore, 
modulo u B there are only two quasienergies, which we designate f2i,2- Finally, we can see 
that this is also the set of energies of the system. Transforming back to the original time 
independent (scalar) gauge, we have 
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-ieEx°t\„,(a) 



[ U W(f)o t (jfe) + u^(t)b\k)] \0)e- iQl 



2t 



in(kd—u> gt) 



u^(t)4 + M W(t)6tl| )e-^' 2 ', 



(A6) 



as the basis set of solutions to the dynamical Schrodinger equation in that gauge (u a,b (t) 
are the parts of the Bloch periodic function u(t) belonging to each of the two bands, in 
an obvious notation). But these must be linear combinations <f>j exp(— iEjt) of energy 
eigenvalue solutions, and we see from (|A6|) that those energies are Ej = Qi } 2 + tiub, for 
arbitary integers n — two interleaved Stark ladders of rung spacing ljb, offset from each 
other by Q = Q 2 — 
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